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Abstract 



Comparing the result of inserting a complete set of physical states in a time 
ordered product of b decay currents with the operator product expansion 
gives a class of zero recoil sum rules. They sum over physical states with 
excitation energies less than A, where A is much greater than the QCD scale 
and much less than the heavy charm and bottom quark masses. These sum 
rules have been used to derive an upper bound on the zero recoil limit of the 
B D* form-factor, and on the matrix element of the kinetic energy operator 
between B meson states. Perturbative corrections to the sum rules of order 
as{A) A^/m^ have previously been computed. We calculate the corrections 
of order as(A) and ^^(A) Pq keeping all orders in A/rricfi, and show that these 
perturbative QCD corrections suppressed by powers of A/nic^b significantly 
weaken the upper bound on the zero recoil B D* form-factor, and also on 
the kinetic energy operator's matrix element. 
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Over the last six years, dramatic progress has been achieved in our understanding of 
exclusive and inclusive B decays. For exclusive decays this resulted from applying heavy 
quark symmetry |I| to relate B decay form-factors and obtain their normalization at zero 
recoil. For example, the form-factors that occur in B ^ D e Ue and B D* e z/g semileptonic 
decays are related by heavy quark symmetry to a single universal function of v ■ v' {v is the 
four- velocity of the B, and v' is that of the recoiling D^*^), and furthermore, this function is 
normalized to unity at zero recoil 

Progress in the theory of inclusive B decays has come from applying the operator product 
expansion and heavy quark effective theory [^] to perform a l/rrib expansion of the time 
ordered product of b decay currents 0. It was found that at leading order in this expansion, 
the inclusive semileptonic B decay rate is equal to the perturbative b quark decay rate. 
There are no nonperturbative corrections at order l/rrif,, and the corrections of order I /ml 
are characterized by only two matrix elements (we use the standard relativistic normalization 
for the B meson states) 

Xi = ^{B{v)\U:HzDfh^:^\B{v)), (1) 

and 

A2 = ^ (Biv) I h^:^ I a,. G^^ h^^ I B{v)) , (2) 

where h^/f> is the b quark field in the heavy quark effective theory • The matrix element 
A2 is scale dependent [l^, and it is determined from the measured B* — B mass splitting, 
A2(mfe) ~ 0.12 GeVl 

Sum rules have been derived that relate exclusive decay form-factors to the matrix ele- 
ments Ai 2 [|n|]. The zero recoil sum rules follow from analysis of the time ordered product 

T,. = ^ Jd'xe-'''-{B\T{Jl{x),JM}\B), (3) 

where J,^ is a 6 — c axial or vector current, the B states are at rest, q = and g° = mh—mc—e. 
Viewed as a function of complex e, T^^, has two cuts along the real e-axis. One, for e ^ 0, 
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FIG. 1. The integration contour C in the complex e plane. The cuts extend to Ree — > ±00. 

corresponds to physical states with a charm quark and the other, for e ^ — 2mc, corresponds 
to physical intermediate states with two h quarks and a c quark. The first cut arises from 
inserting the states between the two currents in the product J, and the second cut arises 
from inserting the states between the currents in the other time ordering J.P . So we arrive 
at 
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The sum over X includes the usual phase space factors, i.e., d^p/2E for each particle in the 
state X. 

Consider integration of the product of a weight function PVaIc) with T^y{e) along the 
contour C shown in Fig. 1. Assuming W is analytic in the shaded region enclosed by this 
contour and averaging over /i = z/ = l,2,3, we get 
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The maximum X mass on the right-hand side of eq. (|^) is determined by where the contour 
C pinches the real axis. For convenience this mass is chosen to be less than 2m;, + to 
prevent the occurrence of states X with b, b, and c quarks. We take the maximum X mass 
to be 2m B- Hereafter it is understood that sums over X only go over states up to mass 
2m B- 

We require that: (i) the weight function Wa be positive semidefinite along the cut so 
that every term in the sum over X on the right-hand side of eq. is non-negative; (ii) 
Wa{0) = 1; (iii) Wa be flat near e = {i.e., at least dWA(e)/de _^ = 0); (iv) and that it 
falls off rapidly to zero for e > A. We want to take A ^ m^b- Then states X other than 
the D* give a contribution to the right-hand side of eq. (|D that is suppressed by (l/mc,fe)^. 
However, in our numerical results we consider A as large as 2 GeV. Although our analysis 
holds for any weight function that satisfies these four properties, for explicit calculations we 
use 

A 2n 

with n = 2, 3, . . . (for n = 1 the integral over e is dominated by contributions from states 
with mass of order m^). These weight functions have poles at e = therefore, as 

long as n is not too large and A is much larger than the QCD scale, Aqcd? the contour in 
Fig. 1 is far from the cut until e is near 2m b. Then we should be able to calculate the integral 
in eq. using the operator product expansion to evaluate the time ordered product. 

The choice of the set of weight functions in eq. is motivated by the fact that for values 

(n) 

of n of order unity all poles of lie at a distance of order A away from the physical cut. 
In this case the integral along the contour C can be computed only assuming local duality 
at the scale 2mB. The dependence of our results on this assumption is extremely weak. 



because for A <^ the weight function is very small where the contour C touches the cut. 
As n — > oo, approaches ^^(A — e) for positive e, which corresponds to summing over all 

hadronic resonances up to excitation energy A with equal weight. Then the poles of 
approach the cut, and the contour C is forced to lie within distance of order A/n from the 



cut at e = A. In this case the evaluation of the integral along the contour C relies also on 
local duality at the scale A.[] 

Neglecting perturbative QCD corrections and nonperturbative effects corresponding to 
operators of dimension greater than five, the operator product expansion gives [H 



1 rjnAA _ _1 _^ (Ai + 3A2)(mb - 3mc) _ 4A2mfc - (Ai + 3A2)(mb - rUc - e) 



(7) 



when J^ = Af, = 07^75 b, and 
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when Jfj, = Vfj^ = C'-ffj^b. Performing the contour integration yields 

77^ E WA[{mx-m,)-{mB-mh)]{27T)^S^{px) {X\A,,\B) 



(8) 
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6^ Sm^mh, 

These equations hold for any that satisfies the four properties mentioned above. Higher 
order terms in the operator product expansion for give contributions with more factors 
of 1/e on the right-hand sides of eqs. (0) and (H). Therefore, if the weight function has 
nonvanishing m'th derivative at e = 0, there are corrections to the right-hand side of eq. 
of order 
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We require that A be large enough compared with the QCD scale Aqcd, so that such terms 
are smaller than those we kept in eq. (Pa|). For m > 1 A can still be smaller than m^fc. Higher 



*In fact, for any sequence of functions analytic in some neighbourhood of the positive real axis 
that converges to ^(A — e), some singularity will approach e = A. Thus, the pinching of the contour 
is inevitable if one uses a weight function that varies rapidly. 
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order terms in the operator product expansion of give corrections to the right-hand side 
of eq. ( pBD of order (Aqcd/"^c,6)^ (Aqcd/^)™^^- This is why we imposed condition (iii). 



For the weight function W^\e) in eq. (|^) the first nonvanishing derivative is at m = 2n. 

We have considered the nonperturbative corrections to the sum rules (^) characterized 
by Ai and There are also perturbative corrections suppressed by powers of the strong 
coupling. These are most easily calculated not in the operator product expansion, but by 
directly considering the sum over states in (|) and replacing the hadronic states by quark 
and gluon states. The perturbative corrections are of two types. There are corrections of 
order a<((mc^b) not suppressed by powers of A/m^fe. These arise, at the parton level, from 
the final state X = c and change the term 1 on the right-hand side of (pa| ) to t]\, where 
riA_ is the usual factor that relates the axial current in the full theory of QCD to the axial 
current in the heavy quark effective theory (at zero recoil). r]A has been calculated to order 
as 0, and terms of order a'j. Po, where Po = (H — I"-/); also known ||l3| , p!^ . Explicitly, 



VA = l In — +o)-—f^0 7-;[ In — + 77 ' W 

TT \mb — rric mb oj tt^ 24 \mb — rric rUb 15/ 

where is the MS coupling evaluated at the scale y/rnt/m^. 

There is another class of perturbative QCD corrections coming from final states X that 
contain a charm quark plus additional partons, e.g., eg, cqq, etc. They give a contribu- 
tion to the right-hand side of equations (§) that is of order [as(A) + . . .] F{A), where the 
ellipses denote terms of higher order in the strong coupling constant as, and for small A, 
F(A) ~ A^/m^j,. We have evaluated the strong coupling constant at the scale A, because 
that characterizes the typical hadronic mass in the sum over X. Note that although these 
corrections are suppressed by powers of A/rric^b they can be as important as the other per- 
turbative corrections we considered, since the strong coupling constant is evaluated at a 
lower scale A. The value for these corrections depends on the precise form of the weight 
function, and we use the ones given in eq. (^). Such perturbative corrections were calculated 
at order as{A) A^/m^ ^ in the limit when the weight function approaches the step-function 



6{A — e) (corresponding to W^'') |jTT],|T5|. As we have already pointed out, the use of such 



b 




SiJlSUUUUL 



b 



FIG. 2. Feynman diagrams that contribute to the order as(A) corrections to the sum rules. 
The black square indicates insertion of the b ^ c axial or vector current. 
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FIG. 3. Feynman diagrams that determine the order a^(A) Pq corrections to the sum rules. 

a weight function relies on local duality at the scale A, so the corrections are expected to be 
less than those stemming from with small n relying on local duality only at the scale 
2mB. We calculated for n > 2 the terms of order as{A) coming from the Feynman diagrams 
in Fig. 2, and the order a^(A) Pq terms arising from the diagrams shown in Fig. 3. Then 
eqs. (01) and (OH) become 



1 



- {ms - m,)] (27r)^ 6'{px) {X\A,\B) 



2 -^2 
Va 2 

mi 



Ai + 3A2 



mi 



mf 



1 



A2 Ai + 3A2 

mi 



Smcmf, 



- (mj, - mi)] {27c)^ SHpx) {X\V,\B) 



TT 



(12a) 



(12b) 



Sm^mf, 



+ f^£(^X<t(A) + Hl(^/J.y«(A). 



TT 



On the right-hand sides of eqs. (|12]) terms suppressed by more than two powers of 
^qcd/"^c,6 or as have been neglected. We have also neglected in ( |12a| ) terms suppressed 
by (Aqcd/"^c,6) (Aqcd/A)^" and in ([T2E| ) terms suppressed by (Aqcd/"^c,6)^ (Aqcd/A)^"-^ 
Perturbative corrections to the terms proportional to Ai,2 are also neglected, and we evalu- 
ate A2 in eqs. (|I^) at the scale m^ (a calculation of QCD corrections to its coefficient would 
resolve this scale ambiguity). For A <^ m^b the functions X*^") and Y^"^^ are given by 



/ ^ /!(")/! 1 



3 Vm^ ml Smcrrib 
6 ml Smcmb 



6 \ml ml SmcmbJ 15 Vm^ ml 3mcmh 



D \m^ mi SmcmhJ 
where the coefficients A^""^ and -B*-"-* are (n > 2) 

= — ^— , = ( + - - In 2V (14) 



nsin(7r/n) ' V2ntan(7r/n) 3 

For A near 1 GeV higher powers of A/m^^fy are important. The analytic expressions for 
X^Z^ and X!^^ are (for A < 2mb) 



-r^(oo) _ ^ , , -.-c./ y-u . --t,/ J p , ^...^ ..-g 
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18ml (A + mc)2 9m; 
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A (A + 2mc) [2(A + rrj-c)^ — 2m^ — (m^ — mc)^] 3m^ — 2mbmc — ml A + 



18m^ (A + mc)^ 9m^ 
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(15) 



{X^ was also calculated in Ref. |TH|. Our result seems to disagree with theirs.) In Figures 



4a and 4b we plot X^°°^ and y versus A using the values mb = 4.8 GeV and mc = 1.4 GeV. 
The thick solid lines are X and the thick dashed lines are Y, while the thin lines are the 
corresponding functions at order A^/m^^. Note that expanding in A/m^b is not a good 
approximation unless A ^ 400 MeV. 

The evaluation of the order al (3q corrections is made relatively simple by the relation 
between the n/ dependent part of the order a] contribution and the order as contribution 
with a finite gluon mass Such a relation holds in the so-called l^-scheme, but throughout 
this paper we present all results in the usual MS scheme. Knowledge of the order al (3q 
corrections allows us to obtain the BLM scale that results from absorbing vacuum 
polarization effects into the running coupling constant. It is generally believed that this 
choice of scale yields a reasonable perturbative expansion. This is also the reason for using 
as{A) in the sum rules in eqs. (p^ . Had we chosen some very different scale /i, the coefficients 
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FIG. 4. X^°°\A) and Y^°°\A) for the a) axial, and b) vector coefficients. Thick sohd fines are 
X wfiile thick dashed fines are Y. The thin sofid and dashed fines are X and Y to order A^/m^^,. 

would contain large logarithms of /i^/A^. Using = 4.8 GeV, = 1.4 GeV, and 
n — s> oo we obtain /Iblm — 0.12 A and fi^hu — 0.17A for the BLM scales for the axial and 
vector current sum rules, respectively. If we demand ^^(/iBLM) < 1, then A needs to be 
above about 3 — 4 GeV, which would completely eliminate the restrictive power of the sum 
rules. Another possibility, which we adopt, is just to use our results as estimates of the 
corrections, but to retain A ~ 1 GeV. Then we find that the order ^^(A) corrections to the 
sum rules are comparable to the order ^^(A) terms, and we can only hope that terms of 
higher order in as(A) are not similarly important. 

The parameters Ai and A2 also occur in the inclusive differential B decay rate, which 
can be expressed in terms of the B and D meson masses and the parameters Ai, A2 and A, 
where 



rriB = mb + A 



Ai + 3A2 



mo = rric + A 



Ai + 3A2 



(16) 



2mf, 2mc 

The pole mass is not a physical quantity, and the perturbative expression for the MS mass 
rrTbirrib) in terms of the pole mass irtb is not Borel summable, giving rise to what is sometimes 



called a "renormalon ambiguity" in the pole mass |18|. However, when the differential 



semileptonic decay rate is expressed in terms of the hadron masses and A, the perturbative 
QCD corrections to the decay rate are also not Borel summable. If A (or equivalently 
the b quark pole mass) extracted from the differential semileptonic decay rate is used to 
get the MS mass these ambiguities cancel, so one can arrive at a meaningful prediction 
for the MS b quark mass. It is fine to introduce unphysical quantities like A as long as 
one works consistently to a given order of QCD perturbation theory and the expansion 
in inverse powers of the heavy quark masses. Since the final results one considers always 
involve relations between physically measurable quantities, any "renormalon ambiguities" 
arising from the bad behavior of the QCD perturbation series at large orders will cancel out 
19| , |20[| . As the left-hand sides of the sum rules in eqs. (0) are physical quantities, the right- 



hand sides, when calculated to all orders in ag, should be free of renormalon ambiguities. 
We checked that the order Aq^j-i/m^ ;, renormalon ambiguity in the quantity ri\ cancels 



against that in the perturbative corrections suppressed by A^/m^f, (such a cancellation was 
conjectured in 



Eqs. (pal) and (pq) have been used to bound the B D* zero recoil form-factor |22,11 



The sum rule ( p.2a| ) implies a bound on the zero recoil B ^ D* matrix element of the axial 
current F^^^,, defined by {D*\Ai\B) = 2yJmD*TTiB Fb~*d' £i, that reads 



Fi^^.<^l--i+^— ^ -, + -2 + 



+ e:A^x^:l(A) + ^p„YSiA). (17) 



Here we used that the contributions of states X of higher mass than the D* to the left-hand 
side of ( |12a| ) are positive, and neglected the very small deviation of W^\{mD* — rric) — 
{niB — rrib)] from unity implied by eq. (||), eq. (piBD , and the relation rriD* — fn^ = 2X2 /nic. 
The positivity of the sum over states X in eq. (|12b| ) implies that 



< 4 - (^) + ^4"^(A) + ^,9„y#(A) . (18) 

This inequality gives a constraint on the heavy quark effective theory matrix element Ai, 
which is strongest when one takes the rric ^ rrih ^ A limit, giving 
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Ai < -3A2 + ^^A2^— + ^^/5oA2^— . 19 

TT 3 TT^ 3 



Neglecting the perturbative corrections suppressed by powers of as(A), eq. ( [I9| ) yields 
Ai < -0.36 GeV^ which in turn implies using eq. (0) that F^^d* < 0.93 (With 
rrib = 4.8 GeV and rric = 1.4 GeV we find rjA = 0.96 following from eq. ([TTl).) To indicate 
the importance of the perturbative corrections proportional to as(A) and aKA) Pq, we give 
the bounds that result when they are included for n = 2, n = 3, and n — oo in Table I. The 
effects of these corrections are smaller if we choose A small (corresponding to suppressing 
the contribution of higher excited states) or if we choose n large (using local duality at the 
scale A). Note that while it is plausible that n can be chosen arbitrarily large as local duality 
is expected to hold at scales much above Aqcd, the relation A ^ Aqcd niust be maintained 
and so A cannot be chosen to be less than about 1 GeV. Using = 3, A = 1 GeV and 
as{l GeV) = 0.45 we obtain the bounds given in Table I. The large magnitude of the second 
order corrections to the sum rules indicates that the series of perturbative corrections might 
be under control only for A significantly above 1 GeV. Such a value for A would greatly 
weaken the restrictive power of the sum rules. Similar comments and conclusions apply to 
two analogous sum rules derived for B* D*^*) transitions in Ref. ||23[ . 



In conclusion, we investigated perturbative corrections to the zero recoil inclusive B 
decay sum rules derived in Ref. |Tl| . We calculated the corrections suppressed by powers of 
A/nic^h at order ^^(A) and order ^^(A) Po corresponding to a set of possible weight functions 





n = 2 


n = 3 


71 = OO 


Ai 


to order Og 


-0.06 GeV^ 


-0.13 GeV2 


-0.17 GeV2 




to order Pq 


0.13 GeV2 


0.06 GeV2 


0.01 GeV2 


Fb^d* 


to order 


0.96 


0.96 


0.95 




to order /5o 


0.99 


0.99 


0.98 



TABLE I. Upper limits on Ai and Fb^d* that can be obtained from the sum rules in eqs. ([T7|) 
and ( [l9| ) with A = 1 GeV. n labels the weight function lU^"^ 
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that determine the contributions of excited hadronic intermediate states. These corrections 
significantly weaken the constraints stemming from the sum rules. It is widely believed that 
Ai < (although in our opinion it has not been proven in QCD for Ai defined by the MS 
subtraction scheme), and we are not aware of any claim that Fb^d* is significantly above 1. 
Due to the size of the A-dependent terms in the sum rules, it is hard to deduce any useful 
model independent bounds. An upper bound below 1 on the zero recoil B D* form-factor 
barely survives these perturbative corrections, and a limit on Ai that restricts it to negative 
values does not. However, it is important to remember that the results in Table I rely on 
the applicability of QCD perturbation theory at a scale A = 1 GeV, and furthermore are 
very sensitive to the value of as at this scale. In the future Ai may be determined from 
experimental data on inclusive B decays [^, and then a bound on Fb^d* that does not 
rely on eq. ([T9| ) can be derived from eq. (p!7|). 

In light of our discussion, we see no reason to think that the original estimates of Fb^d* , 
based on model calculations, the structure of terms arising at order l/m^f^ [^, and on 
chiral perturbation theory [Q, badly underestimated the corrections. Our results 

show that the zero recoil sum rules do not demand a larger deviation of Fb^d* from rjA, 
even if the D* does not saturate the sum over states X. We cannot prove at this point that 
such a deviation does not occur. However, in the absence of any such indication, it is most 
natural to think that Fb^d* — Va = 0.96 holds to an accuracy of about the canonical size 
of the corrections, that is within (500 MeV/2mc)^ — 3%. 
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